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64c.  Uvw  triple  products  presented  in  wall-coordinates  at  Station  3. 
64d.  Uvw  triple  products  presented  in  wall-coordinates  at  Station  4. 
64e.  Uvw  triple  products  presented  in  wall-coordinates  at  Station  5. 

*  t  * 

tfLLL 


101 

101 

102 

103 

103 

104 

104 

105 

105 

106 
106 

107 

108 
108 
109 

109 

110 
110 
111 
111 
112 
113 

113 

114 

114 

115 

115 

116 
116 

117 

118 
118 
119 

119 

120 


Figure  64f.  uvw  triple  products  presented  in  wall-coordinates  at  Station  6 .  120 

Figure  64g.  uvw  triple  products  presented  in  wall-coordinates  at  Station  7 .  121 

Figure  64h.  uvw  triple  products  presented  in  wall-coordinates  at  Station  8 .  121 

Figure  64i.  uvw  triple  products  presented  in  wall-coordinates  at  Station  9 .  122 

Figure  65a.  IF  triple  products  presented  in  wall-coordinates  at  Station  1 .  123 

Figure  65b.  IF  triple  products  presented  in  wall-coordinates  at  Station  2 .  123 

Figure  65c.  IF  triple  products  presented  in  wall-coordinates  at  Station  3 .  124 

Figure  65d.  IF  triple  products  presented  in  wall-coordinates  at  Station  4 .  124 

Figure  65e.  IF  triple  products  presented  in  wall-coordinates  at  Station  5 .  125 

Figure  65f.  I F  triple  products  presented  in  wall-coordinates  at  Station  6 .  125 

Figure  65g.  IF  triple  products  presented  in  wall-coordinates  at  Station  7 .  126 

Figure  65h.  IF  triple  products  presented  in  wall-coordinates  at  Station  8 .  126 

Figure  65i.  IF  triple  products  presented  in  wall-coordinates  at  Station  9 .  127 

Figure  66a.  IF  triple  products  presented  in  wall-coordinates  at  Station  1 .  128 

Figure  66b.  IF  triple  products  presented  in  wall-coordinates  at  Station  2 .  128 

Figure  66c.  IF  triple  products  presented  in  wall-coordinates  at  Station  3 .  129 

Figure  66d.  IF  triple  products  presented  in  wall-coordinates  at  Station  4 .  129 

Figure  66e.  IF  triple  products  presented  in  wall-coordinates  at  Station  5 .  130 

Figure  66f.  IF  triple  products  presented  in  wall-coordinates  at  Station  6 .  130 

Figure  66g.  IF  triple  products  presented  in  wall-coordinates  at  Station  7 .  131 

Figure  66h.  F  triple  products  presented  in  wall-coordinates  at  Station  8 .  131 

Figure  66i.  IF  triple  products  presented  in  wall-coordinates  at  Station  9 .  132 

Figure  67a.  IF  triple  products  presented  in  wall-coordinates  at  Station  1 .  133 

Figure  67b.  IF  triple  products  presented  in  wall-coordinates  at  Station  2 .  133 

Figure  67c.  IF  triple  products  presented  in  wall-coordinates  at  Station  3 .  134 

Figure  67d.  IF  triple  products  presented  in  wall-coordinates  at  Station  4 .  134 

Figure  67e.  IF  triple  products  presented  in  wall-coordinates  at  Station  5 .  135 

Figure  67f.  IF  triple  products  presented  in  wall-coordinates  at  Station  6 .  135 

Figure  67g.  IF  triple  products  presented  in  wall-coordinates  at  Station  7 .  136 

Figure  67h.  IF  triple  products  presented  in  wall-coordinates  at  Station  8 .  136 

Figure  67i.  IF  triple  products  presented  in  wall-coordinates  at  Station  9 .  137 
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ABSTRACT 


Experimental  data  obtained  in  a  two-dimensional  turbulent  boundary  layer  (TBL) 
flow  and  a  three-dimensional  (3D)  pressure-driven  TBL  flow  around  a  wing/body  junction 
are  presented.  Reynolds  number  effects  are  investigated.  Fine  spatial  resolution  laser- 
Doppler  velocimeter  measurements  of  mean  velocity  and  higher  order  statistics  including 
the  third  order  products  measured  for  9  stations  are  reported.  Reynolds  number  based 
on  momentum  thickness  (f?e6)  of  the  developing  2-D  boundary  layer  on  the  wall  was 
Re6 =23200,  and  the  approach  nominal  reference  velocity  of  the  undisturbed  airflow  was 
32  m/sec.  The  wing  shape  used  had  a  3:2  elliptical  nose,  NACA  0020  profiled  tail  and  it 
was  mounted  on  a  flat  wall.  Distribution  of  measured  static  pressure  variation  on  the  wall, 
turbulent  kinetic  energy  and  some  derived  quantities  such  as  the  flow  angle,  flow  gradient 

angle,  shear  stress  angle,  and  modeling  parameters  such  as  x/v2, 
v+yt+vw2),  B2  =177 [( u v2)2 +( v2  w)2]  are  also  discussed. 


NOMENCLATURE 


U/UjU, 

U,V,W 

Ue 

8 

cP 

P 

P* 

Pstag 

Pref 

y 

y+ 

p 

V 

U2,'v2,'w2 

uv,Vw,Vw 


Reynolds’  stress  tensor 

Triple  product  tensor 

mean  velocity  components  (m/sec) 

mean  velocity  magnitude  above  the  layer  edge  at  each  station 
Boundary  layer  thickness.  Distance  where  the  velocity  magnitude  in 
the  layer  reaches  to  0.995  of  the  Ue. 

pressure  coefficients 

mean  static  pressure  measured  on  the  wall-surface 
pressure  fluctuating  component 
stagnation  pressure 

reference  static  pressure  of  the  undisturbed  freestream  measured  at 
the  tunnel  entrance, 
distance  from  wall  (microns) 

nondimensional  distance  to  the  wall,  y+=yujv 
flow  density  (kg/m3) 
kinematic  viscosity  (m2/sec) 

Reynolds’  normal  stresses  (m2/sec2) 

Reynolds’  shear  stresses  (m2/sec2) 

Skin-friction  ve\oc\ty=JxJ p 
Wall  shear  stress  (Newton/m2) 


A 

TKE 

B 

b2 

Pw 

Pfs 


shear  stress  magnitude=\/(i/v)2+(  vw)2 

_  ^(ZTv)2+(vw)2 


Townsend’s  structural  parameter,  At=. 


u‘  +  v*  +  w 


Turbulent  kinetic  energy=.l  (u2+v2+'w2) 


B=  v7/ ( JFv+v2 + vw2) ,  triple  product  ratio 

B2=  v7/[(uv2)2+(v2w)2],  triple  product  ratio 
wall-stress  angle  with  respect  to  the  tunnel  coordinates 
free-stream  angle  with  respect  to  the  tunnel  coordinates 


CHAPTER  I.  INTRODUCTION 


Wing-body  junction  flows  have  been  the  subject  of  many  studies  in  the  past.  The 
flow  at  the  nose  of  the  wing  has  aperiodic  (bimodal)  double-peaked  fluctuating  velocity 
histograms  with  very  large  Reynolds  stress  values  compared  to  the  approaching  boundary 
layer  values  (Devenport  and  Simpson,  1990a).  Shinpaugh  and  Simpson  (1995)  used  a 
scanning  LDV  and  a  pressure  transducer  to  further  investigate  this  velocity  and  the 
pressure  field  simultaneously. 

The  approach  boundary  layer  separates  due  to  the  adverse  pressure  gradient  in 
this  region  and  rolls  in  towards  the  floor,  generating  backflow.  The  flow  visualization  water 
tunnel  studies  of  Khan  et  al.  (1995)  and  Kim  et  al.,  (1991)  also  show  these  characteristics 
and  describe  the  time  development  of  the  motion.  The  vortical  flow  is  stretched  around 
the  wing  and  is  reoriented  by  the  wing,  which  results  in  vortical  flow  along  the  sides  of  the 
wing-body  junction. 

The  3-D  TBL  of  the  pressure-driven  type  has  been  the  subject  of  many  papers. 
Reviews  of  recent  experimental  data  were  done  by  Anderson  and  Eaton  (1 989),  Johnston 
and  Flack  (1996),  Schwarz  and  Bradshaw  (1994),  and  Olgmen  and  Simpson  (1992, 
1993).  DNS  studies  of  Spalart  (1989)  and  Moin  et  al.  (1990)  also  show  the  capability  of 
CFD  in  explaining  the  flow  physics  at  low  Reynolds  numbers.  Three-dimensional 
boundary  layer/vortex  interaction,  which  is  relevant  to  the  work  undertaken  here,  was 
studied  by  Eaton’s  group  (Eaton,  1995). 

In  the  present  study  the  detailed  mean  flow,  Reynolds  stress  and  triple  product 
characteristics  of  3-D  turbulent  boundary  layers  forming  on  the  wall  around  a  wing/body 
junction  flow  at  9  stations  are  presented.  The  wing  shape  used  had  a  NACA  0020  tail;  3:2 
elliptical  nose  wing  (maximum  thickness  t=  7.17  cm).  Figure  1  shows  the  measurement 
locations.  All  the  data,  unless  otherwise  specified,  are  presented  in  tunnel  coordinates 
where  x  is  along  the  tunnel  centerline  and  y  axis  is  perpendicular  to  the  wall.  Table  1 
gives  some  of  the  flow  parameters  at  the  measurement  locations. 

Flow  around  this  wing  at  a  lower  approach  Reynolds  number  of  f?e0= 5940  is 
extremely  well  documented  and  detailed  references  can  be  found  in  the  work  by 
Simpson  (1995, 1996)  and  Olgmen  and  Simpson  (1995a,  1996e).  Me  Mahon,  Merati  and 
Yoo  (1987),  Me  Mahon,  Hubbart  and  Kubendran  (1982),  and  Dickinson  (1986)  made 
measurements  with  a  similar  model  at  other  Reynolds  numbers  to  report  the  flow  field 
around  and  in  the  wake  of  the  wing  using  hot-wire  probes. 

The  present  study  differs  from  these  previous  studies  with  more  detailed  and  more 
accurate  near  wall  information  of  not  only  of  the  mean  flow  and  Reynolds  stresses  but 
also  of  the  triple  products.  Present  detailed  data  are  the  highest  Re6  wing/body  junction 
flow  studied.  In  high  Re6  flows  there  is  a  dire  need  for  data  sets  which  include  triple 
products  for  better  understanding  of  the  diffusion  processes,  and  data  sets  should  also 
include  near  wall  sublayer  data  since  the  most  interesting  flow  physics  occur  in  that 
region. 

The  five-component  laser-Doppler  velocimeter  that  was  used  for  the  near-wall 
measurements  and  the  data  reduction  techniques  were  previously  discussed  (Olgmen 
and  Simpson,  1995a,  b).  The  effective  probe  volume  had  a  diameter  of  30  microns.  The 
probe  volume  could  traverse  as  close  as  50  microns  to  the  wall  and  the  wall  could  be 


found  with  an  uncertainty  of  ±10  microns.  In  this  study  the  five  component  LDV  system 
was  modified  to  reach  higher  distances  away  from  the  wall.  The  probe  volume  could  be 
traversed  to  6  inches  away  from  the  wall  with  reduced  spatial  resolution.  The  near  wall 
measurements  were  made  using  the  original  design,  and  measurements  above  40mm 
were  made  using  the  modified  LDV  system.  Two  overlapping  measurements  made  at 
each  station  were  used  to  form  composite  profiles. 

In  Chapter  II,  design  modifications  of  the  5  component  LDV  system  are  described. 
The  high  Ree  flow  generation  procedure  is  described  in  Chapter  III.  The  present  2D  zero- 
pressure  gradient  high  Ree  data  are  compared  to  previous  high  Ree  2DTBL  data.  In 
Chapter  IV  selected  3-D  data  are  presented  and  the  differences  in  flow  characteristics  at 
selected  locations  due  to  different  flow  conditions  are  explained.  Unless  otherwise 
specified  the  data  are  presented  in  tunnel  coordinates.  The  report  includes  the  mean 
velocity,  Reynolds’  stress,  all  triple  products.  The  derived  quantities  such  as:  skewnesses 

for  different  velocity  components,  turbulent  kinetic  energy,  u vfu1,  shear  stress  magnitude, 

Townsend’s  structural  parameter  A1f  (JF+wz)fvz,  and  velocity  magnitude  are  presented 
in  Chapter  V.  In  Chapter  VI  the  high  Reynolds  number  data  are  compared  to  previously 
obtained  low  Ree=5940  data  in  the  same  tunnel  using  the  same  wing  geometry  and  with 
a  reference  nominal  velocity  of  Uref=27.5  m/sec. 


CHAPTER  II. DESIGN  MODIFICATIONS  FOR  THE 
5  COMPONENT  LDV  SYSTEM 


In  the  first  part  of  this  Chapter  the  design  modifications  applied  to  the  5  component 
LDV  system  are  described,  and  in  the  last  part  the  data  reduction  program  is 
summarized.  The  design  modifications  were  required  to  traverse  the  measurement  probe 
volume  in  a  range  0  to  6.35“  away  from  the  wind  tunnel  floor.  Expected  boundary  layer 
thicknesses  were  less  than  6.35".  The  design  changes  included  changing  the 
transmitting/receiving  optical  lenses,  and  changing  the  orientation  of  the 
transmitting/receiving  optical  heads  (see  figure  2).  The  focal  lengths  of  the 
transmitting/receiving  optics  were  changed  since  the  probe  volumes  were  required  to  form 
further  away  from  the  optics  heads.  The  orientation  of  the  transmitting  optics  heads  were 
changed  to  allow  the  traversing  without  major  changes  to  the  system. 

The  modified  system  has  56  degrees  between  the  bisectors  of  the  beams 
emerging  from  the  two  transmitting  probe  heads,  as  opposed  to  90  degrees  obtained  with 
the  original  LDV  system  (figure  2).  This  was  accomplished  by  inserting  17°  apex  angled 
aluminum  wedges  between  the  transmitting  optics  heads  and  the  translators  that  hold  the 
heads  in  their  place.  In  the  modified  system  the  velocities  are  thus  not  measured  in  an 
orthogonal  frame,  however  the  three  components  measured  could  easily  be  transformed 
into  an  orthogonal  frame. 

The  transmitting  lenses  used  in  the  modified  system  had  200mm  (JML  Direct  Inc., 
part  no:  DBL14130)  focal  lengths.  The  separation  of  the  beams  emerging  from  the 
transmitting  head  was  kept  as  9.1mm,  which  resulted  in  an  included  angle  of  2.606° 
between  the  beams  passing  through  the  focusing  lenses.  With  a  nominal  effective  focal 
length  of  200  mm,  a  1  mm  diameter  laser  beam  focuses  down  to  174.7pm  in  diameter. 
The  length  of  the  crossing  along  the  optical  axis  is  calculated  as  7.681  mm,  and  the  fringe 
spacing  for  the  514.5nm  laser  beams  is  calculated  as  1 1 .3  pm.  The  light  collection  probe 
volume  for  the  system  is  defined  by  the  overlapping  portions  of  the  individual  probe 
volumes.  Thus  the  long  interference  region  of  individual  probe  volumes  do  not  reduce  the 
spatial  resolution.  However  the  probe  volume  of  the  system  in  the  modified  system  is 
larger  than  the  original  5CLDV  due  to  the  larger  size  of  the  focused  beam  diameters.  On 
the  other  hand  the  modified  system  is  intended  to  be  used  only  in  regions  of  the  flow 
away  from  the  wall  where  the  velocity  gradients  are  not  large,  therefore  the  reduced 
spatial  resolution  does  not  effect  the  overall  quality  of  the  measured  velocities. 

The  receiving  lens  unit  uses  two  achromatic  lenses  with  175  mm  (JML  Direct  Inc., 
part  no:  DBL14110)  and  60  mm  (JML  Direct  Inc.,  part  no:  DBL13985)  focal  lengths, 
instead  of  two  120mm  focal  length  achromats  used  in  the  original  system.  The  1:1 
magnification  ratio  of  the  previous  system  is  thus  reduced  to  2.92:1  magnification  ratio. 
The  receiving  fiber  core  diameter  was  kept  as  50  microns.  Previous  LDV  measurements 
with  the  original  5CLDV  system  showed  that  the  light  was  collected  only  through  the  30 
pm  portion  of  the  fiber  core.  Based  on  this  fact  with  a  2.92:1  magnification  ratio  the  light 
is  collected  from  88  pm  portion  of  the  overlapping  measurement  probe  volumes.  The 
receiving  fiber  thus  collects  the  light  from  7-8  fringes. 

Besides  the  mechanical  differences  between  the  original  and  the  modified  systems, 
specific  electronic  data  acquisition  parameters  needed  to  be  adjusted  for  the  modified 


system.  In  the  modified  design  Macrodyne  3100  frequency  domain  processors  used 
needed  64  record  length  setting/frequency  analysis.  The  record  length  is  the  number  of 
samples  used  by  the  processors  to  calculate  the  Doppler  frequency  using  fast-Fourier 
transform  process.  The  setting  proved  to  be  important,  since  the  data  obtained  with  32 
record  length  with  a  larger  than  95%  validation  rate  was  observed  to  result  in  very  high, 
unusual  normal  stress  values.  The  record  length  setting  did  not  effect  the  shear  stress 
information  obtained  form  the  data.  The  noise  content  in  separate  processor  do  not 
correlate  with  each  other  thus  the  shear  stresses  are  not  affected.  For  the  original  system 
the  record  length  was  kept  as  32. 

The  measurements  with  the  near  wall  LDV  probe  were  repeated  up  to  6  times  at 
different  stations.  The  outer  layer  profiles  were  also  acquired  many  times.  The  main 
reason  was  to  explain  the  differences  observed  between  the  measurements  obtained  with 
the  near  wall  and  outer  layer  LDV  probes.  The  discrepancy  in  the  measured  velocity 
values  with  original  and  the  modified  systems  were  minimized  by  using  the  correct  record 
length  setting  for  the  Macrodyne  frequency  domain  processors.  The  uncertainties 
presented  are  calculated  using  differences  in  the  measured  quantities  from  different  data 
sets  obtained  at  the  same  station  (Table  2.). 

II  A.  Summary  of  data  reduction  program 

In  this  section  the  data  reduction  program  is  summarized.  The  detailed  data 
reduction  program  description  can  be  found  in  the  paper  by  Olgmen  and  Simpson, 
(1 995b).  At  every  measurement  point  total  number  of  30000  samples  with  more  than  95% 
validation  ratio  of  the  Macrodyne  units  were  acquired.  Each  of  the  three  simultaneous 
components  measured  were  used  to  form  three  histograms  in  the  optics  coordinate 
system.  The  histograms  were  first  cleaned  in  that  coordinate  system  by  fitting  parabolas 
to  the  logarithm  of  the  histogram  values  one  on  each  side  of  the  peak  and  then  discarding 
the  points  which  laid  outside  the  parabolas  crossing  with  the  ln(1)  line.  Once  a  component 
of  the  velocity  was  discarded  the  accompanying  two  velocity  measurements  were  also 
discarded.  Next  the  measured  simultaneous  velocity  components  were  used  to  calculate 
the  simultaneous  velocity  components  in  tunnel  coordinates  by  a  transformation.  The 
cleaning  process  was  repeated  one  more  time  in  the  tunnel  coordinates  to  eliminate  the 
noise  at  the  skirts  of  the  histograms.  The  noise-free  data  were  used  to  calculate  new 
histograms  using  the  inverse-velocity-magnitude  as  the  weighting  function.  The  bias- 
corrected  histograms  were  next  used  to  calculate  the  mean  velocity  and  higher  order 
terms.  Table  2  gives  21:1  uncertainties  calculated  using  the  Kline  and  Me  Clintock 
method.  The  uncertainties  in  individual  terms  were  estimated  using  2  to  6  separate  data 
sets  measured  at  each  measurement  location. 


CHAPTER  III.GENERATION  OF  THE  Re,= 23200 
APPROACH  FLOW  AND  COMPARISON  WITH  DATA 

FROM  LITERATURE 


111  A.  Roughness  plate  and  tunnel  flow  characteristics  comparison 

Two-dimensional  roughness  elements  were  used  to  artificially  thicken  the  boundary 
layer.  The  configuration  of  these  elements  (0.125"  square  ribs)  is  shown  in  the  figure  3a. 
The  leading  edge  of  this  configuration  is  placed  at  a  distance  of  53"  from  the  entrance  to 
the  test  section.  The  test  section  is  over  29  ft  long  (about  353  inches)  with  a  width  of  3 
•  feet  and  a  variable  height.  The  test  section  is  divided  into  three  sections  numbered  from 

the  contraction  as  sections  one,  two,  and  three.  The  first  section  has  an  accelerating  flow 
up  to  five  feet  from  the  contraction  and  approximately  a  constant  inviscid  core  velocity 
downstream.  The  second  section  has  a  constant  height,  but  the  third  section  has  a 
varying  height  that  has  been  suitably  adjusted  to  get  a  zero  pressure  gradient  2-D 
^  boundary  layer.  The  third  section  begins  at  a  distance  of  about  200.375  inches  from  the 

entrance  to  the  test  section.  Plexiglass  inserts  of  0.25"  thick  were  placed  along  the  side 
walls  beginning  at  the  streamwise  distance  of  200.37"  to  a  distance  of  260.31".  The 
heights  (H)  of  this  section  at  various  streamwise  distances  (X)  are  shown  in  the  following 
table.  The  dimensions  are  in  inches 


X(in) 

200.38 

206 

212 

218 

224 

230 

236 

242 

248 

H(in) 

10.875 

10.85 

10.875 

10.906 

10.938 

10.938 

10.938 

10.938 

10.969 

X(in) 

254 

260 

266 

272 

278 

284 

290 

296 

308 

H(in) 

11.031 

11.031 

10.938 

10.844 

10.812 

10.75 

10.688 

10.75 

10.75 

e  rest  of  the  test  section  downstream  is  at  a  constant  height  of  1 0.75  inches. 


height  adjustment  in  this  third  section  was  done  to  obtain  a  nominally  zero  pressure 
boundary  layer.  The  free  stream  velocity  in  the  test  section  as  a  function  of  the 
streamwise  distance  (X)  is  shown  in  figure  3b.  The  high  Reynolds  number  measurements 
were  made  at  the  streamwise  location  of  276  inches,  approximately.  Appendix  I  describes 
in  detail  the  additional  pitot-static  tube  spanwise  measurements  made  in  boundary  layers 
to  document  the  symmetry  of  the  flow.  Measurements  were  also  used  to  validate  the 
boundary  layer  flow  quality. 

HI  B.  DNW  data  (Nockemann,  et.  al,  1994  )  and  the  present  2D  data  comparison 

DNW  (f?e6= 22140,  and  Ree  =20920)  and  the  present  2D  data  (f?ee  =23200)  mean 
velocity  and  the  normal  stresses  are  plotted  in  figures  4  and  5,  respectively.  The  mean 
U /ux  values  differ  =2.5  from  each  other  throughout  the  measurement  range.  The 
differences  between  the  nondimensional  normal  stresses  are  approximately  0.5.  The 
differences  are  not  easily  attributable  to  the  size  of  the  hot-wire  probe  or  the  technique 
itself.  The  single  normal-wire  probe  used  at  DNW  had  length/diameter=220  with  a 
corresponding  l+=22.5  (for  Ree  =20920  case,  l+=47  for  the  Ree=41260  case).  These 


values  are  close  to  the  values  suggested  by  Ligrani  and  Bradshaw  (1987)  l/d>200  and 
l+=<20  for  turbulence  statistics  to  be  free  of  wire  length  effects.  The  DNW  skin-friction 
values  were  measured  using  a  Preston  tube,  and  a  semi-empirical  relationship  by 
Femholz  (with  a  ±8%  variation  in  the  calculated  skin-friction  coefficient  values  due  to 
different  methods  used  (Fernholz  et  al.,  1995).  For  the  present  data,  using  the  semi-log 
law-of-the-wall  concept  with  Coles’  constants  (0.41 ,  and  5)  resulted  in  ux= 1 .08  while  with 
Femholz’s  constants  (0.4,  and  5.1)  the  ut=1.057.  For  sublayer  data  using  a  function  fit 

in  the  form  of  U=a(y+c)-b(y+c)4 ,  where  a=u2/v,  ux=  0.98  was  calculated.  The  skin  friction 
velocity  was  chosen  as  the  average  of  these  numbers  as  uT=1.03.  The  ux= 1 .03  used  in 
the  present  study  is  lower  than  the  1 .057  calculated  using  Fernholz’s  constants.  Thus  U+ 
values  are  higher  than  the  law-of-the-wall  variation.  The  small  size  of  the  probe  volume 
used  in  these  studies  eliminates  the  probe  volume  size  effects  in  the  turbulence 
measurements.  The  30|j.m  near  wall  probe  diameter,  and  88|im  outer  layer  probe 
diameter  correspond  to  l+=1.85  and  l+=5.43,  respectively,  well  below  the  requirements 
defined  by  Ligrani  and  Bradshaw  (1987). 


Ligrani  and  Bradshaw  (1987)  show  that  with  a  probe  size  1  <  l+  <  10  the 


changes  semi-logarithmically  from  a  value  2.72  at  ReB= 1000  to  3.03  at  Ree =10000.  For 
higher  Re6  flows  hot  wire  probes  can  not  be  made  small  enough  for  conventional  size 
tunnels,  which  restricts  the  use  of  the  technique  at  high  Ree  flows.  The  maximum  value 
of  the  ratio  drops  at  higher  Ree  flows  with  the  increased  value  of  l+.  However,  the  present 


=3.1 3  calculated  at  Red =23200  is  along  the  semi-logarithmic  variation  observed 


by  Ligrani  and  Bradshaw  (1987. 


The  mean  velocity  and  the  —  normal  stress  profiles  for  DNW,  Re  =20920  case 


ur 


and  the  present  2D  data  follow  each  other  closely  (figures  12  and  13).  The  agreement 

IF 

between  the  —  normal  stresses  is  better  above  y+=200  as  compared  to  Ree=22140 

Ur 


case.  The  DNW,  Ree  =20920  and  the  present  flow  U  mean  velocity  values  are  also  in 
better  agreement  than  the  flee  =22140  case.  The  reason  for  this  is  probably  a  combined 
effect  of  the  hot-wire  probe  size  and  the  skin-friction  velocity  estimation  methods  used 
both  for  the  DNW  and  the  present  study. 

The  DNW  data  obtained  at  Ree  =20920  (Femholz  and  Finley,  1996,  their  figure  47) 


match  well  above  y+=2000  for  the  —  (figure  5).  Their  data  vary  semi-logarithmically  from 

Ur 

=1.2  at  y+=30  to  =1.45  at  y+=1500,  while  the  present  data  quickly  rise  to  =1.6  at  y+=100 
and  increases  slightly  till  y+=1250.  The  peak  value  location  at  y+=1250  is  closely  predicted 
by  the  relation  defined  by  Fernholz  and  Finley  (1996)  for  the  maximum  location  as 


U/6 

The  —  peak  values  presented  by  Fernholz  and  Finley  are  *2.5  (their  fig.  51) 

Ur 

-rp. 

while  the  peak  value  for  the  present  data  is  *2.9  (figure  5).  For  the  present  case  — 

Ur 

varies  semi-logarithmically  from  y+=100  to  1500,  similar  to  the  variation  observed  by 
Fernholz  and  Finley  for  f?ee=21390  case. 

The  DNW  Uv  data  show  a  peak  around  *-1.0  at  y+=400  (figure  6).  The  variation 
above  y+=1000  is  about  the  same  as  the  present  data. 

The  DNW  triple  product  data  and  the  present  triple  product  data  are  plotted  in 

figures  7  and  8.  The  low  -  u  observed  at  y+=250  for  the  present  data  is  not  observed 

Ur 

by  Fernholz  and  Finley.  Bruns  Ree= 2573  data  presented  by  Fernholz  and  Finley  show 
a  dip  at  this  y+,  however  the  values  are  not  less  than  zero.  Similarly  large  values  of 

Ur 

observed  for  the  present  data  at  y+=300  is  not  observed  for  data  presented  by  Fernholz 
and  Finley. 

The  "u  skewness"  near  wall  variation  presented  by  Fernholz  and  Finley  for 
Ree  =20920  case  is  very  similar  to  the  variation  of  the  present  data  till  y+=30  (their  fig.  66), 
however  the  positive  values  of  "u  skewness"  and  following  semi-logarithmic  variation  are 
not  observed  (figure  9).  Data  presented  by  Fernholz  and  Finley  (1996)  show  that  the  "u 
skewness"  near  the  wall  in  the  sublayer  is  independent  of  ReQ.  Above  y+=1000  the 
present  data  and  the  DNW  data  for  Ree= 22140  follow  each  other  closely.  The  large 
variation  at  the  layer  edge  is  due  to  the  presence  of  turbulent-nonturbulent  interface. 

The  “v  skewness"  values  presented  here  are  similar  to  the  data  presented  by 
Fernholz  and  Finley  (their  fig.  72),  however  the  near  wall  negative  values  observed  for 
the  present  data  are  not  observed  for  the  other  data  sets  (figure  10).  The  "sweep"  and 
"ejection"  type  of  motions  observed  near  the  wall  requires  that  the  "v  skewness"  to  be 
less  than  zero  around  the  y+=15  where  the  turbulent  kinetic  energy  production  is 
maximum.  Thus  it  is  believed  that  the  present  data  is  more  correct  indicator  of  the  "v 
skewness"  variation  near  the  wall. 

The  "w  skewness"  values  for  both  flows  are  close  to  zero  (figure  11). 

Figures  presented  in  this  section  clearly  show  that  the  DNW  turbulence  data  are 
much  more  scattered  than  the  present  data.  The  DNW  data  uncertainties  are  also  not 
given  in  the  papers  documenting  the  data.  In  conclusion  the  present  data  set  is  a  data 
set  with  low  uncertainties  (Table  2)  and  much  less  scattered  than  the  DNW  data  set. 


CHAPTER  IV.EXPERIMENTAL  RESULTS  FOR  THE 
HIGH  REYNOLDS  NUMBER  WING/BODY  JUNCTION 

FLOW 


IV  A.  Oil  flow  visualization  pictures 

Oil  flow  visualization  pictures  obtained  in  two  different  Re%  number  flows  are  presented 
in  figure  12a  and  12b.  The  figures  show  that  the  incoming  flow  separates  from  the  wall, 
rolls  inward  towards  the  wing/wall  junction  at  the  nose  region  and  follows  the  contour  of 
the  wing  to  form  a  horse-shoe  shape.  Near  the  leading  edge  between  the  separation  line 
and  the  leading  edge  there  is  another  line  where  the  limiting  streamlines  pass  through. 
This  line  is  formed  due  to  the  chaotic  nature  of  the  flow  field  and  it  is  named  as  line-of- 
low-shear,  since  the  shear  stresses  following  this  line  is  smaller  than  the  surrounding 
regions.  The  flow  separates  near  the  trailing  edge  to  form  the  fish-tail  separation  region. 
The  figure  13  shows  the  schematic  of  the  oil  flow  visualization  and  some  important 
distances  which  were  measured  from  the  oil-flow  pictures.  These  lengths  are  given  in  the 
following  table. 


Reynolds 

number 

V 

;xsep 

X|ois 

D 

L 

W 

W 

¥*max 

5940 

-0.45 

-0.27 

1.38 

3.76 

1.72 

2.02 

23200 

-0.475 

0.201 

3.77 

1.84 

1.994 

The  abbreviations  "sep"  stands  for  separation,  nlol 

feir  stands  for  line  of  low  shear. 

IV  B.  Mean  static  surface-pressure  distribution 

The  wall  static  pressure  on  the  tunnel  floor  was  measured  using  a  plexiglass  plate 
with  230  pressure  taps  around  the  contour  of  the  wing  and  a  small  metal  plate  with  139 
pressure  taps  at  the  nose  region.  A  Scanivalve  system  (CTLR2P/S2-S6)  was  connected 
to  all  taps.  Two  Setra  239  pressure  transducers  were  used  to  measure  the  port  static 
pressure  and  reference  dynamic  pressure.  A  data  acquisition  board  (DT2801)  in 
conjunction  with  a  PC-A7286  was  used  to  acquire  the  data  from  the  pressure  transducers. 
Once  a  port  was  connected  to  the  Scanivalve,  pressure  measurement  was  taken  after  a 
settling  period  of  15  seconds. 

The  figure  14  is  a  contour  plot  of  the  Cp  values.  The  wing’s  presence  generates 
longitudinal  adverse  pressure  gradient  region  around  the  nose  region  of  the  wing.  The 
pressure  gradients  become  favorable  following  the  contour  of  the  wing  towards  the 
maximum  thickness  location.  Near  the  trailing  edge  the  adverse  pressure  gradients  are 
generated.  Details  of  surface  mean-static  pressure  measurements  are  given  in  Appendix 
II. 


IV  C.  Mean  Flowfield 


The  mean  velocity  profiles  are  presented  in  figures  15  and  16.  The  U  mean 
velocity  show  logarithmic  variation  regions  similar  to  2-D  flow  U  mean  velocity  variations. 
In  2-D  flow  the  logarithmic  regions  for  different  Reynolds  number  flows  overlap  up  to 
y+=400  to  form  a  "law-of-the-wall"  region.  However  in  a  3-D  flow  the  profiles  at  different 
stations  do  not  overlap.  There  does  not  exist  a  similarity  law-of-the-wall  for  3-D  TBL  flows 
(Olgmen  and  Simpson,  1992).  The  U  mean  velocity  first  decelerates  until  station  4  and 
then  accelerates. 

The  W  mean  velocity  profiles  show  that  the  peak  values  in  the  profiles  are  located 
further  away  from  the  wall  proceeding  downstream.  The  pressure  gradients  are  most 
effective  on  the  low  momentum  lateral  flow  nearer  the  wall  and  the  turning  of  the  flow  is 
propagated  throughout  the  layer.  In  the  present  flow  the  lateral  pressure  gradient  is 
positive  at  upstream  stations,  negative  at  stations  4-7,  and  then  positive  again.  Since 

(d2Wldy2)  is  proportional  to  this  spanwise  pressure  gradient,  the  changes  in(9W/9y) 
propagate  away  from  the  wall.  Thus,  the  maximum  magnitude  of  the  profiles  increase  until 
station  5  and  than  decrease  to  form  the  wave-like  velocity  profiles  at  station  9.  The  outer 
layer  W  velocity  gradient  is  practically  zero  indicating  a  zero  flow  gradient  angle,  and 

reduced  productions  for  the  w2,  Uw  and  vw  stresses.  This  also  implies  that  the  presence 
of  the  wall  is  propagated  out  into  the  layer  at  a  much  further  downstream  distance  than 
a  lower  Reynolds  number  flow. 


IV  P.  Stress  variation 


The  Reynolds-stress  transport  equation  with  summation-subscript  notation  can  be 
written  as  (Daly,  and  Harlow,  1970): 
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The  —  profiles  at  the  first  two  stations  show  that  the  fluctuations  of  the  flow  are 
Ur 

highly  increased  throughout  the  layer,  especially  in  the  logarithmic  variation  region  of  the 
U  mean  velocity  (figures  17a  and  17b).  This  is  however  followed  by  a  very  high  reduction 
in  the  fluctuations  in  the  outer  layers  at  the  stations  where  the  longitudinal  pressure 
gradient  is  favorable  (starting  with  station  4).  The  dimensional  data  show  that  adverse 

longitudinal  pressure  gradient  reduces  the  IP  value  below  y+=100,  and  increases  the 
values  above  this  y+  location.  The  pressure  gradient  affects  the  flow  near  the  wall  and 
above  y+=100  in  the  opposite  sense.  The  near  wall  values  below  y+=100  are  increased 


at  stations  4  through  7  where  the  longitudinal  pressure  gradients  are  negative,  the  flow 

is  accelerating,  and  production  of  TP  increases.  However,  above  y+=100  the  values  are 
reduced  at  the  same  stations,  suggesting  that  there  are  2  distinct  regions  within  the 

layers.  In  this  region  the  production  of  the  TP  is  much  lower  than  nearest  the  wall  and  the 
rapid  distortion  takes  place. 

TP 

The  —  values  are  also  increased  compared  to  the  2-D  values  at  the  first  2 

'A 

stations  where  the  longitudinal  pressure  gradient  is  an  adverse  pressure  gradient  (figures 
18a  and  18b).  However  the  values  are  smaller  than  the  2-D  values  proceeding 

downstream  until  station  8.  On  the  contrary  the  dimensional  data  show  that  the  TP  normal 
stress  values  are  overall  increased  throughout  the  layers  with  the  three  dimensionality  of 
the  flow  as  compared  to  the  2D  flow  data. 

The  —  stresses  are  also  increased  at  the  first  2  stations  and  they  gradually  relax 

ut2 


until  station  8,  which  may  be  due  to  the  redistribution  of  elevated  TP  values  (figures  19a 

and  19b).  In  a  2-D  flow  the  production  term  for  the  TP  is  practically  zero,  and  the  main 

source  for  the  TP  is  due  to  the  redistribution  of  TP.  The  values  increase  at  stations  8  and 

9.  Similar  to  the  TP  normal  stress,  the  TP  normal  stress  values  increase  proceeding 
downstream,  once  the  data  are  plotted  in  dimensional  units.  All  the  three-dimensional 
values  presented  in  dimensional  units  are  larger  than  the  2D  data  values. 

The  Tv  shear  stress  values  are  increased  in  the  outer  layers  of  the  first  2  stations 
compared  to  2-D  values  and  the  maximum  stress  locations  are  further  away  from  the 
wall.  This  is  similar  to  the  observations  made  by  other  researchers  in  2-D  flows  under 

adverse  pressure  gradients  (figures  20a  and  20b).  However  the  Tv  values  reduce  quickly 
below  2-D  values  with  the  accelerating  flow  in  the  outer  regions.  Below  y+=1 00  the  values 

increase  proceeding  downstream.  The  dimensional  Tv  magnitude  values  are  increased 
with  the  three  dimensionality  of  the  flow  in  the  outer  layers  above  y+=300.  Near  wall 
dimensional  values  however  show  a  dependence  on  the  longitudinal  pressure  gradient 
sign.  With  favorable/adverse  pressure  gradients  the  magnitude  increases/decreases. 

The  production  term  for  the  Tw  stress  is  -(Tv  BW/dy+vw  BU/By)  and  the 

production  term  for  the  Tw  stress  is  -(TP  BWIBy+Tw  BV/By).  The  production  of  these 

stresses,  especially  that  of  Tw,  are  close  to  zero  where  the  dW/dy= 0  (figures  21a  to 
22b).  The  peak  value  locations  and  the  sign  of  the  profiles  follow  the  variation  of  the 
d  W/By.  Both  stresses  show  wave  type  nature  throughout  the  layers.  Even  though  the 

value  of  the  Tw  is  the  largest  it  drops  out  of  the  equations  used  for  computations  once 
the  boundary  layer  approximations  are  applied. 

In  an  attempt  to  relate  the  variation  of  the  U  mean  velocity  and  the  TP  normal 

stress  the  following  relation  was  observed.  Dimensional  data  showed  that  the  TP  normal 
stress  and  the  mean  velocity  are  inversely  proportional  to  each  other.  In  order  to  describe 


12 


this  co-dependent  variation  the  product  (u^XU)  was  calculated  using  the  tunnel 

coordinates  data.  The  (TPX U)  product  nondimensionalized  with  the  Uref3  are  shown  in 
figures  (figures  23a  and  23b).  The  product  profiles  overlap  with  each  other  above  5mm, 
which  roughly  corresponds  to  the  region  above  0.04  *  5  (boundary  layer  thickness).  The 
product  can  be  interpreted  as  the  faster  mean  flow  results  in  smaller  velocity  fluctuations, 
or  faster  flow  reduces  the  longitudinal  normal  stress  values. 

In  another  attempt  to  describe  the  co-dependent  variation  of  the  U  mean  velocity 

and  the  TP  normal  stress  in  tunnel  coordinates,  Batchelor’s  (1 953)  "Sudden  distortion  of 
a  turbulent  stream"  concept  was  used.  The  method  assumes  that  the  distortion  of  the 
stream  takes  place  so  rapidly  that  the  inertia  and  the  viscous  forces  arising  from  the 
turbulent  motion  have  no  effect.  With  these  assumptions,  the  problem  becomes  linear. 
It  is  also  assumed  that  in  the  regions  upstream  and  the  downstream  of  the  distorting 
section  of  the  stream  tube,  the  mean  flow  is  uniform  and  the  turbulence  is  homogenous. 

The  equations  describing  the  co-dependent  variation  of  the  IP  normal  stress  is  than  given 
as: 
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For  the  lateral  components  the  following  equation  holds: 
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and,  the  equation  for  the  TKE  ratio  reads  as: 
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At  all  stations  IP  decreases  in  the  outer  layer  approximately  as  given  by  "p".  The  use  of 
"p"  computed  using  the  ratio  of  the  local  mean  velocity  magnitudes  at  2  streamwise 

locations  for  "c"  also  collapsed  the  TP  data  within  10%  for  y  >  104  microns  (figures  24a 

and  24b).  The  use  of  Batchelor’s  "v"  for  the  distortion  of  the  lateral  ~P  and  TP  stresses 
and  the  TKE  ratio  did  not  collapse  those  data  (figures  25  and  26). 


IV  E.  Triple  Products 


Triple  products  are  presented  in  dimensional  (m/s)3  and  nondimensional  correlation 
coefficient  form  in  tunnel  coordinates  in  figures  27a  to  36b.  Presentation  in  correlation 
coefficient  form  shows  the  combined  effect  of  the  mean  velocity  variations  due  to  the 
pressure  gradients  both  on  the  normal  stresses  and  triple  products. 

The  TPv/  u2v'  (Fig.  27b)  shows  the  transfer  of  the  I P  normal  stress  with  the  “v" 

fluctuations.  The  TP  on  the  average  diffused  towards  the  wall  below  y+=20  and  away  from 
the  wall  above  this  height.  Maximum  TKE  production  occurs  at  y+=15  to  20,  and  most  of 


the  TKE  is  stored  in  TF.  The  dimensional  plot  of  the  quantity  (Fig.  27a)  shows  that  all  the 
station  profiles  collapse  on  top  of  another  above  1 .5  cm  away  from  the  wall  and  the 
variation  is  a  semi-logarithmic  variation.  This  indicates  that  the  turbulent  diffusion  process 
is  not  affected  by  the  three  dimensionality  of  the  flow  in  the  outer  region  of  the  flow  in  this 
high  Reynolds  number  flow.  The  effect  of  three  dimensionality  is  felt  further  from  the  wall 
in  this  high  Re  flow  than  the  previous  study  of  the  same  configuration  in  a  lower  Ree 
(=5940)  flow  because  the  current  boundary  layer  is  much  thicker. 

The  TFw  variations  (Fig.  28a)  are  large  only  near  the  wall  below  y+=100.  This 
shows  that  the  turbulent  structures  are  almost  aligned  in  the  same  direction  of  the  mean 

flow  at  the  first  station,  a  similar  conclusion  observed  from  the  TFv  profiles.  This  is  also 

clearly  observed  in  the  TF  variation  (Fig.  36a).  The  triple  products  involving  "w" 
fluctuations  above  y=1 .5  cm  are  not  affected  by  the  three  dimensionality  of  the  flow;  the 
turbulent  diffusion  process  is  similar  to  2D  flow  diffusion.  Three-dimensional  effects  have 
not  propagated  into  the  outermost  part  of  this  thick  boundary  layer. 

The  variations  of  the  Fw2  (Fig.  31a),  TF  (Fig.  34a),  and  the  IFF  (Fig.  30a)  triple 

products  are  similar  to  each  other.  The  TF,  TF,  and  the  TF  normal  stresses  are 
transferred  on  the  average  in  the  same  fashion  by  the  "u"  fluctuating  velocity  component. 
Below  y+=100  the  transfer  is  associated  with  the  higher-than-the-mean  velocity 

fluctuations  and  above  that  height  the  transfer  is  with  lower  speed  fluctuations.  PositiveTF 

and  negative  TFv  nearer  the  wall  indicates  "sweep"  type  of  motion  (+u,-v)  is  the  dominant 
motion  in  this  region.  Above  y+=100,  the  ejection  type  of  motion  (-u,+v)  is  dominant.  Semi- 

logarithmic  variations  of  the  TF/u' 3  (Fig.  34b)  and  the  TivF/u’w 2  (Fig.  31b)  show  that 
between  y+=1 00-4000  these  triple  products  are  scaled  with  the  normal  stresses.  The 

dimensional  profiles  of  the  TF  and  of  the  IFF  triple  products  overlap  on  top  of  another 

for  each  station,  respectively.  The  nondimensional  TF  (Fig.  35b),  TTF  (Fig.  32b)  plots 

show  that  the  transfer  of  the  TF  and  the  TF  normal  stresses  are  away  from  the  wall 
approximately  above  y+=25. 


CHAPTER  V  DERIVED  QUANTITIES 


V  A.  Flow.  Row  Gradient,  and  Shear  Stress  Angles 

The  flow  angle,  flow  gradient  angle  and  shear  stress  angles  in  tunnel  coordinates 
are  shown  in  figures  37  to  45.  The  near  wall  flow  angle  first  decreases  until  station  3  and 
then  increases  at  the  downstream  stations  to  become  positive  at  station  8.  Previous 
three-dimensional  experimental  data  show  that  the  shear  stress  angle  lags  behind  the 
flow  gradient  angle.  However  at  stations  1  and  2  in  the  high  Reynolds  number  case  the 
shear  stress  angle  near  the  wall  seems  to  lead  the  flow  gradient  angle;  the  shear  stress 
angle  has  a  larger  magnitude  than  the  flow  gradient  angle  at  these  stations.  At 
downstream  stations  the  flow  gradient  angle  leads  the  shear  stress  angle.  The  changes 
in  the  flow  gradient  angle  results  in  gradual  change  of  the  shear  stress  angle  in  the  region 
above  y+=1 00.  Above  y+=1 000  the  shear  stress  angle  seems  to  change  very  slightly  away 
from  the  station  1  mean  free-stream  flow  direction.  The  changes  in  the  mean  flow 
direction  due  to  pressure  gradients  and  the  subsequent  changes  in  the  shear  stress 
values  does  not  seem  to  affect  the  shear  stress  angle.  The  turbulent  structure  direction 
above  y+=1000  does  not  seem  to  change  with  the  three-dimensionality  of  the  flow  in  this 
high  Ree  flow. 


V  B.  Townsend’s  Structural  Parameter.  A, 


The  A,  =\(uv)2  +(vw)2  / (IF  +~v2  +~w2) ,  Townsend’s  structural  parameter  is  shown 
in  figure  46.  The  maximum  value  at  the  first  three  stations  is  0.1  which  is  lower  than  0.14 
observed  for  the  2-D  station.  The  peak  value  obtained  in  the  outer  regions  increase  to 
0.14  proceeding  downstream.  The  data  show  that  the  shear  stress  magnitude  and  the 
turbulent  kinetic  energy  of  the  flow  at  the  first  4  stations  in  the  3-D  flow  are  not  as  well 
correlated  as  in  the  2-D  flow. 

V  C.  Turbulent  Kinetic  Energy 

The  Turbulent  Kinetic  Energy  (TKE={1F +v* +w2)/2)  profiles  are  shown  in  figure 
47.  This  figure  shows  that  the  TKE  of  the  flow  stays  constant  above  y+=1000.  Large 
changes  occur  nearer  the  wall  where  the  production  of  the  TKE  is  maximum. 

V  D.  Stress.  Normal  Stress  Ratio,  t /IF 


The  x/v 2  profiles  show  that  the  shear  stress  magnitude  and  the  IF  normal  stress 
are  highly  correlated  with  each  other  (figure  48).  The  ratio  is  practically  a  constant  around 
0.6  within  the  y+=1 00  to  1 000  range.  This  fact  was  previously  shown  using  9  different  data 
sets  (Olgmen  and  Simpson,  1995). 


V  E.  Triple  Product  Ratio.  B 


The  parameter  B=\F/(u2v+\F  +VW2 )  is  the  ratio  of  the  perpendicular-to-the-wall  velocity 
fluctuation  triple  product  to  the  triple  product  terms  which  appear  in  the  TKE  transport 
equation  (figures  49a  and  49b).  The  parameter  is  defined  to  reduce  the  number  of 
independent  terms  that  need  to  be  modeled  in  governing  equations.  "B"  varies  as  a 
constant  with  a  value  between  0.2  and  0.4  between  y+=300  to  3000.  The  "B"  parameter 
was  also  observed  to  be  within  0.2  to  0.4  using  other  available  3-D  data  sets  (Olgmen  et 
al.  1995). 


V  F.  Triple  Product  Ratio.  B„ 

The  parameter  B2 = Vs/ [  ( u  v2)2  +  ( \F  w)2 ]  is  the  ratio  formed  among  the  triple 

products  which  appear  in  the  transport  equations  for  the  Vs  normal  stress  (figures  50a 
and  50b).  The  ratio  is  approximately  a  constant  (ratio~1)  between  y+=500  to  6000,  similar 
to  the  variations  observed  for  other  3-D  flow  data  (Olgmen,  et  al.,  1998). 

V  G.  Z/v  /  IF.  Shear  Stress  to  Normal  Stress  Ratio 

The  Uv/TF,  shear  stress  to  normal  stress  ratio  which  shows  the  correlation 
between  these  stresses  is  plotted  in  Figure  51 .  The  V  fluctuating  velocity  is  usually 
assumed  to  be  unaffected  by  the  large  longitudinal  inactive  motions  as  described  by 

Bradshaw.  Thus  the  Uv  should  not  be  affected  by  these  inactive  motions,  since  the 
Fourier  frequency  of  the  large  scale  motions  generating  the  inactive  motion  and  the 
Fourier  frequency  of  the  small  scale  structures  generating  the  V  fluctuating  velocity  are 
largely  separated  from  one  another.  The  quantity  therefore  indicates  two  possible  flow 
phenomena:  1)  Either  "u"  and  "v"  fluctuating  velocities  are  less  correlated  at  stations 
where  the  pressure  gradient  is  adverse  or  2)  the  inactive  motion  is  more  pronounced  at 
adverse  pressure  gradient  stations. 

V  H.  (IF  +  W}  !_  yh  Ratio  of  Normal  Stresses 

The  ratio  of  the  stresses,  {TF  +  wF)  /  is  shown  in  figures  52a  and  52b.  For  an 
isotropic  turbulent  flow  the  ratio  is  equal  to  2.  However  the  ratio  varies  semi- 
logarithmically  above  y+=65  in  the  present  flow  (figure  52a).  The  ratio  near  the  layer  is 
less  than  2,  indicating  the  anisotropic  distribution  of  the  turbulent  kinetic  energy  (figure 
52b).  Figure  52b  also  shows  that  the  semi-logarithmic  variation  is  a  first  order 
approximation  to  the  variation  above  y+=65.  The  values  show  a  sharper  decay  in  the  high 
Reynolds  number  case  compared  to  the  f?e9= 5940  case.  The  Ree=5940  case  decays 
from  a  value  of  =4.5  at  y+=100  to  =3.5  at  y+=1000,  compared  to  the  variation  from  =6  to 
=3.5  in  the  same  range  for  the  ReQ= 23200  case. 


V  I.  /  \fzF  .  Ratio  of  Square  Root  of  Normal  Stresses 
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The  ratio  y  v2  /  V ^  is  shown  in  Figure  53.  The  ratio  shows  the  anisotropy  of  the 
TKE  distribution  between  the  longitudinal  and  perpendicular  normal  stresses.  At  adverse 
pressure  gradient  stations  the  ratio  is  lower  than  the  2D  values  indicating  the  increased 
anisotropy  of  the  flow.  However,  the  distributions  follow  a  similar  shape  at  all  stations.  In 
comparison  to  DNW  (German-Dutch  windtunnel)  data  (Fig.  63  of  Femholz  and  Finley, 
1996),  the  2D  station  variation  of  the  parameter  is  very  much  similar  to  their  data 
obtained  for  the  Ree= 20920  case.  The  values  change  semi-logarithmically  from  0.5  to 
0.58  in  the  y+=100  to  1000  range.  The  value  is  increased  to  0.85  by  y+=6000. 


V  J.  \l\flr  /  vtr  ,  Ratio  of  Square  Root  of  Normal  Stresses 


The  profiles  of  the  ratio  yw2  / 


are  shown  in  Figure  54.  Similar  to  the 


/  v Tr  ratio,  the  ywr  / 


i  u*  also  show  the  anisotropy  of  the  normal  stresses.  The 
variation  of  the  2D  station  data  of  the  present  study  is  very  similar  to  the  DNW  data 
measured  at  f?e6  =20920  (Femholz,  Finley,  1996,  Fig.  65).  In  the  y+=100  to  1000  range 
the  ratio  is  approximately  a  constant  around  0.65.  The  values  are  increased  to  0.85  at 
y+=6000. 


CHAPTER  VI  STUDY  OF  REYNOLDS  NUMBER 

EFFECTS 


In  this  Chapter  three-dimensional  data  sets  obtained  in  high  (Ree= 23200)  and  low 
(f?ee=5940)  Reynolds  number  flows  are  compared.  Both  data  sets  were  obtained  in  the 
same  tunnel  using  the  same  wing  geometry.  The  low  Re%  flow  data  have  been  presented 
previously  by  Olgmen,  and  Simpson  (1995a). 

VI  A.  Mean  Flowfield 


The  mean  velocity  profiles  are  presented  in  figures  55.  Due  to  the  similar 
distributions  of  the  wall  static  pressure  the  mean  flow  field  near  the  wall  vary  similarly  in 
both  Reynolds  number  cases.  The  U  mean  velocity  has  a  larger  logarithmic  region 
extending  further  away  from  the  wall  similar  to  variations  observed  in  2-D  high  Reynolds 
number  flows.  The  flow  decelerates  until  station  4  and  then  accelerates. 

The  W  mean  velocity  in  the  first  four  stations  are  skewed  in  one  direction  only.  At 
stations  5  and  downstream,  the  pressure  gradient  sign  changes  and  the  profiles  show 
minima.  The  near  wall  low  momentum  fluid  region  is  the  region  of  the  layers  first  affected 
by  the  pressure  gradients.  This  results  in  negative  velocity  components  away  from  the 
wall  in  the  wall-coordinates  system. 

VI  B.  Stress  variation 


Except  at  stations  3  and  4  the  —  profiles  for  the  high  and  the  low  Reynolds 

Ur 

number  flows  follow  each  other  very  closely  below  y+=1000,  indicating  that  the  wall  skin- 

friction  is  closely  related  to  the  v2  normal  stress  (figure  56  and  57).  At  stations  3  and  4 
the  high  Reynolds  number  values  are  lower  than  the  low  Reynolds  number  profile  values. 

IF 

The  —  values  at  every  station  for  the  high  Reynolds  number  case  is  higher  than 

Ur 

the  low  Reynolds  number  profiles.  Near  the  wall  this  might  be  an  indicator  of  higher 
inactive  motion  for  the  higher  Reynolds  number  case.  The  maxima  location  in  the  layers 
both  for  the  high  and  the  low  Reynolds  numbers  cases  are  at  the  same  y+  value  away 

from  the  wall.  The  production  term  for  the  TF  normal  stress  is  -(2  Uv—) .  The  uv 

dy 

profiles  nondimensionalized  using  the  skin-friction  velocity  differ  from  each  other  away 
from  the  wall  while  the  nondimensional  U  mean  velocity  have  the  same  slope  for  both  the 

Reynolds  number  cases.  The  higher  level  of  production  may  be  a  reason  for  the  elevated  IF 
values  in  the  outer  region  of  the  profiles.  Another  reason  maybe  due  to  the  larger 
boundary  layer  thickness.  The  diffusion  process  would  take  a  longer  distance  farther 
downstream  for  its  effects  to  be  felt  in  the  outer  layer.  The  viscous  diffusion  and  the 
viscous  dissipation  terms  are  large  only  near  the  wall. 


The  values  of  —  are  also  larger  throughout  the  layers  including  station  7,  which 

u r2 

may  be  due  to  the  redistribution  of  elevated  U1  values.  In  a  2-D  flow  the  production  term 
for  the  TF  is  practically  zero,  and  the  main  source  for  the  TF  is  due  to  the  redistribution 

of  TF.  For  the  high  Reynolds  number  case  even  in  the  2-D  flow  the  T F/uT2  is  larger,  which 
supports  this  view. 

The  variations  of  the  shear  stress  profiles  are  similar  to  each  other  in  both 
Reynolds  number  cases,  especially  below  y+=100.  The  nondimensional  stress  values  are 
larger  for  higher  Reynolds  number  case  above  y+=100  as  expected,  since  the  boundary 

layers  extend  to  higher  y+  units.  The  nondimensional  Uw  and  vw  profiles  follow  each 
other  more  closely  until  y+=1000  than  the  Uv  stress.  At  station  6  the  maximum  value  of 
the  Uw/u T2  is  larger  in  magnitude  than  the  maximum  value  of  the  uv/tif.  However,  the 

importance  of  the  Uw  is  diminished  once  the  boundary  layer  approximations  are  applied 
to  the  transport  equations  for  the  Reynolds’  stresses,  since  this  stress  drops  out  of  the 
equations  to  be  modeled. 


VI  C.  Triple  products 

Triple  products  at  selected  stations  are  presented  in  the  correlation  form  (figures 
58  through  67).  The  triple  product  correlation  in  wall  variables  at  both  Reynolds  numbers 
follow  each  other  within  the  uncertainty  bands  nearer  the  wall,  y+  <  1000.  Obviously  the 
triple  product  values  change  with  different  three-dimensional  flow  conditions,  and  the 
values  depend  on  the  coordinate  system  studied.  However,  for  these  two  flows  the 
Reynolds  number  does  not  seem  to  play  a  role  in  the  development  of  the  triple  products 
for  30  <  y+  <  1000.  The  triple  product  values  for  the  higher  Re  case  extend  to  large  y+ 
further  out  in  the  layers  since  the  layers  are  much  thicker. 


CHAPTER  VII.  CONCLUSIONS 


A  three-dimensional  high  Reynolds  number  flow  around  a  wing/body  junction  has 
been  experimentally  investigated  using  laser-Doppler  velocimetry.  The  mean  velocities, 
shear  stresses,  and  10  triple  products  as  well  as  the  pressure  distribution  on  the  wall  are 
presented.  Some  important  conclusions  are  as  follows:  The  triple  products 
nondimensionalized  using  the  normal  stresses  follow  each  other  at  all  stations  above 
y+=100.  The  shear-stress  angle  lags  behind  the  flow-gradient  angle  at  station  3  and 
downstream  stations.  Parameters  B  and  B2  are  constants  in  the  outer  layer  also  for  this 
high  Reynolds  number  flow.  Townsend’s  structural  parameter  A1  is  not  constant 
throughout  the  layers. 

Two  separate  high  Reynolds  number  data  sets  (Ree=5940  and  Ree  =23200) 
around  a  wing/body  junction  are  compared  with  each  other.  For  the  U  mean  velocity  uz 
is  shown  to  be  the  scaling  velocity.  The  \N/ux  mean  velocity  variations  although  close  to 
each  other  show  that  the  magnitude  of  the  velocity  for  the  Ree  =23200  case  is  smaller. 

The  TF/u?,  and  w^/u?  normal  stresses  are  higher  for  the  Ree =23200  case  especially  in 

the  outer  regions  of  the  layers.  The  normal  stress  does  not  seem  to  depend  on 
Reynolds  number.  The  shear  stresses  normalized  with  the  skin-friction  velocity  show  that 
there  are  two  different  regions  forming  in  this  type  of  flow:  A)  below  y+=100  the  stresses 
tend  to  overlap  with  each  other,  B)  the  magnitude  of  the  normalized  shear  stresses  are 
higher  for  the  high  ReQ  case,  in  the  outer  layers.  Triple  product  profiles  for  30  <  y+  <  1000 

seem  to  be  independent  of  the  Ree .  Future  publications  will  discuss  the  physical  meaning 
of  these  data  in  more  detail. 
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TABLE  1.  LASER-DOPPLER  VELOCIMETER  LOCATIONS  AND  FLOW  PARAMETERS. 
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-0.2643 


-0.2544 


-0.1599 


-0.1094 


Maximum  thickness  of  the  wing  (t)  =  7.17  cm,  Nominal  reference  velocity  (Ure()  =  32m/sec, 

Reynolds  number  based  on  momentum  thickness  in  2D  flow  without  the  wing’s  presence  (Ree)  =  23200. 
Wing  chord  length  =  30.48  cm. 


TABLE  2.  21:1  ODDS  ±2(7  UNCERTAINTIES  OF  MEAN 
VELOCITIES,  REYNOLDS’  STRESSES  AND  TRIPLE 
PRODUCTS. 
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TABLE  3.  SOME  LENGTH  SCALES  OBTAINED 
IN  TUNNEL  COORDINATES. 
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Ue=Velocity  magnitude  at  the  layer  edge. 
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Figure  la.  Wing  shape  and  measurement  locations.  Full  arrows  wall-stress  direction. 
Empty  arrows  free-stream  direction. 


Figure  1  b.  Schematic  of  wing  and  measurement  locations. 
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Figure  2.  Short  and  long  range  transmitting  lens  configuration  for  the  5CLDV.  The 
receiving  lens  configuration  is  for  the  long  range  design. 
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Freestream  Velocity  Vs 
Streamwise  Distance 
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Figure  3b.  Free-stream  velocity  variation  along  the  tunnel  axis. 
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Figure  4.  DNW  (F?e0=2214O,  f?e0  =20920)  data  and  present  2D  (Re0  =23200)  data  mean 
velocity  comparison. 
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Figure  5.  DNW  (/?e0  =22140,  ReQ  =20920)  data  and  present  2D  (f?e0  =23200)  data  normal 
stress  comparison. 


0 


3i 


n/Azn 


Figure  7.  DNW  (Ree  =22140)  data  and  present  2D  (Red  =23200)  data  triple  product 
comparison. 
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Figure  8.  DNW  (Ree  =22140)  data  and  present  2D  (f?ee  =23200)  data  triple  product 
comparison. 
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Figure  10.  DNW  (/?e6= 22140)  data  and  present  2D  ( Ree  =23200)  data  "v  skewness' 
comparison. 


Figure  10a.  DNW  (f?ee= 22140)  data  and  present  2D  (f?e9  =23200)  data  “v  skewness" 
comparison. 


3 


1  j  ■  mr»TH  I  I  I  I  ll|  •  i  1 1  mi  I  II  III  j  IHIH'I  I  I  I  I  llj  ■  I  ■  Mill  I  I  I  I  III  'III 


,~J  [  ■  '  1 1  "I  ■  I  I  I  1 1 1 1  . 1  t  I  1 1 1 1 1  .  .  ■  I...I  I  I  1 1 1 1 1  •  ■  •  ■ . 

100  2  3  101  2  3  102  2  3  1 03  2  3  10+  2 

y* 


Figure  11.  DNW  (Ree= 22140)  data  and  present  2D  (Ree  =23200)  data  "w  skewness' 
comparison. 
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Figure  11a.  DNW  {Ree  =22140)  data  and  present  2D  {Red  =23200)  data  "w  skewness1 
comparison. 


Figure  12a.  Oil  flow  visualization  picture  for  Ree= 5940  case 
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Figure  12b.  Oil  flow  visualization  picture  for  Re6  =23200  case. 
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Figure  14.  Wall  static  pressure  coefficient  distribution. 
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Figure  16a.  W,  dimensional  mean  velocity  profiles  in  tunnel  coordinates 
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Figure  16b.  W/uT  mean  velocity  profiles  in  tunnel  coordinates 
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Figure  17b.  IF/142  normal  stress  profiles  in  tunnel  coordinates 
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Figure  18b.  v^/i/T2  normal  stress  profiles  in  tunnel  coordinates 


Figure  19a.  TP,  dimensional  normal  stress  profiles  in  tunnel  coordinates. 


Figure  21b.  uw/u T2  shear  stress  profiles  in  tunnel  coordinates. 
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Figure  23a.  Nondimensional  U  X  IP  (mean  velocity  X  normal  stress)  vs.  "y"  distance 
away  from  the  wall  for  the  high  Reynolds  number  (fte6 =23200)  case. 
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Figure  23b.  Nondimensional  UX?  (mean  velocity  X  normal  stress)  vs.  “y"  distance 
away  from  the  wall  for  the  low  Reynolds  number  (/?ee= 5940)  case. 
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Figure  24a.  TF  normal  stress  nondimensionalized  by  Batchelor’s  V  parameter. 


Figure  24b.  The  IF  normal  stress  nondimensionalized  by  Batchelor’s  "jj."  parameter. 
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Figure  26.  The  IF+v^+w2,  sum  of  normal  stresses  nondimensionalized  by  Batchelor’s 
"tke  multiplier"  parameter. 
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Figure  27a.  Dimensional  triple  product  in  tunnel  coordinates.  Symbols  denote  the 
measurement  stations. 
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Figure  27b.  Triple  product  in  tunnel  coordinates.  Symbols  denote  the  measurement 
stations. 
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Figure  28a.  Dimensional  triple  product  in  tunnel  coordinates.  Symbols  denote  the 
measurement  stations. 
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Figure  28b.  Triple  product  in  tunnel  coordinates.  Symbols  denote  the  measurement 
stations. 


Figure  29a.  Dimensional  triple  product  in  tunnel  coordinates.  Symbols  denote  the 
measurement  stations. 
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Figure  29b.  Triple  product  in  tunnel  coordinates.  Symbols  denote  the  measurement 
stations. 
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Figure  30a.  Dimensional  triple  product  in  tunnel  coordinates.  Symbols  denote  the 
measurement  stations. 


Figure  30b.  Triple  product  in  tunnel  coordinates.  Symbols  denote  the  measurement 
stations. 


Figure  31a.  Dimensional  triple  product  in  tunnel  coordinates.  Symbols  denote  the 
measurement  stations. 


Figure  31b.  Triple  product  in  tunnel  coordinates.  Symbols  denote  the  measurement 
stations. 


Figure  32b.  Triple  product  in  tunnel  coordinates.  Symbols  denote  the  measurement 
stations. 
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Figure  33b.  Triple  product  in  tunnel  coordinates.  Symbols  denote  the  measurement 
stations. 
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Figure  34a.  Dimensional  triple  product  in  tunnel  coordinates.  Symbols  denote  the 
measurement  stations. 
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Figure  34b.  Triple  product  in  tunnel  coordinates.  Symbols  denote  the  measurement 
stations. 
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Figure  35a.  Dimensional  triple  product  in  tunnel  coordinates.  Symbols  denote  the 
measurement  stations. 
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Figure  35b.  Triple  product  in  tunnel  coordinates.  Symbols  denote  the  measurement 
stations. 
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Figure  36a.  Dimensional  triple  product  in  tunnel  coordinates.  Symbols  denote  the 
measurement  stations. 


Figure  36b.  Triple  product  in  tunnel  coordinates.  Symbols  denote  the  measurement 
stations. 
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re  37.  Flow  angle,  flow-gradient  angle,  and  shear-stress  angle  calculated  at  Station 


Figure  38.  Flow  angle,  flow-gradient  angle,  and  shear-stress  angle  calculated  at  Station 
2. 


station  3 


Figure  39.  Flow  angle,  flow-gradient  angle,  and  shear-stress  angle  calculated  at  Station 
3. 
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Figure  40.  Flow  angle,  flow-gradient  angle,  and  shear-stress  angle  calculated  at  Station 
4. 
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Figure  41 .  Flow  angle,  flow-gradient  angle,  and  shear-stress  angle  calculated  at  Station 
5. 


Figure  42.  Flow  angle,  flow-gradient  angle,  and  shear-stress  angle  calculated  at  Station 
6. 
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Figure  44.  Flow  angle,  flow-gradient  angle,  and  shear-stress  angle  calculated  at  Station 
8. 
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FA,  FGA,  SSA 


Figure  46.  Townsend’s  A1  parameter  for  high  Re  flow. 
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Figure  48a.  1/S=x/vs  flow  parameter. 
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Figure  48b.  VS=x/v2  flow  parameter  calculated  for  different  Reynolds  number  data 
obtained  at  Station  5. 
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Figure  48c.  1/S=x/v^  flow  parameter  calculated  for  different  Reynolds  number  data 
obtained  at  Station  5. 
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B2=v3/  ((v2u)2+(v2w)2)1/2 


Figure  50a.  B2 =vsl[(uvzf +( Vs w)2\  flow  parameter. 
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Figure  50b.  B2  =  v*![ [uv1)2  +( v*  w)2]  flow  parameter. 
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Figure  52a. 
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Figure  52b.  — _ — ,  ratio  of  normal  stresses. 

v2 


Figure  53.  v’/u’,  ratio  of  the  rms  of  the  fluctuating  velocities. 


Figure  54.  w’/u\  ratio  of  the  rms  of  the  fluctuating  velocities. 
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e  55a.  U  and  W  mean  velocities  presented  in  wall-coordinates  at  Station  1 . 
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Figure  55b.  U  and  W  mean  velocities  presented  in  wall-coordinates  at  Station  2. 
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Figure  55c.  U  and  W  mean  velocities  presented  in  wall-coordinates  at  Station  3. 
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Figure  55d.  U  and  W  mean  velocities  presented  in  wall-coordinates  at  Station  4. 
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Figure  55g.  U  and  W  mean  velocities  presented  in  wall-coordinates  at  Station  7. 
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Figure  55h.  U  and  W  mean  velocities  presented  in  wall-coordinates  at  Station  8. 
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Figure  55i.  U  and  W  mean  velocities  presented  in  wall-coordinates  at  Station  9. 
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Figure  56b.  Normal  stresses  presented  in  wall-coordinates  at  Station  2. 
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•  Figure  56d.  Normal  stresses  presented  in  wall-coordinates  at  Station  4. 


79 


80 


81 


U2/ Ur  2,V2/ Ur  2,  W2/ U 


Figure  57a.  Reynolds’  stresses  presented  in  wall-coordinates  at  Station  1 . 


Figure  57b.  Reynolds’  stresses  presented  in  wall-coordinates  at  Station  2. 
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Figure  57d.  Reynolds’  stresses  presented  in  wall-coordinates  at  Station  4. 
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Figure  57e.  Reynolds’  stresses  presented  in  wall-coordinates  at  Station  5. 
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Figure  57f.  Reynolds’  stresses  presented  in  wall-coordinates  at  Station  6. 
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e  57g.  Reynolds’  stresses  presented  in  wall-coordinates  at  Station  7. 
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Figure  57h.  Reynolds’  stresses  presented  in  wall-coordinates  at  Station  8. 


Figure  57i.  Reynolds’  stresses  presented  in  wall-coordinates  at  Station  9. 
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Figure  59i.  TFw  triple  products  presented  in  wall-coordinates  at  Station  9. 
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Figure  60e. 
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Figure  61  d. 
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Figure  61f. 
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Figure  65b. 
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Figure  65d. 


Figure  65f. 
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roducts  presented  in  wall-coordinates  at  Station  2. 
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Figure  66c. 
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Figure  66i.  v5  triple  products  presented  in  wall-coordinates  at  Station  9. 
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Figure  67f. 
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Figure  67i.  w5  triple  products  presented  in  wall-coordinates  at  Station  9. 


Appendix  I 

'Pitot  Probe  Boundary  Layer  Measurements" 
by  Cindy  Whitfield 


introduction 

Pressure  measurements  were  made  across  a  turbulent  boundary  layer  using  a  pitot  probe 
in  the  Virginia  Tech  Boundary  Layer  Wind  Tunnel.  Velocity  profiles  were  calculated  from 
the  acquired  pressure  data  for  11  spanwise  locations  in  the  tunnel  test  section.  The 
momentum  thickness  (0),  and  the  corresponding  Ree  were  calculated  from  the  velocity 
profiles.  These  results  verified  that  the  installed  roughness  plate  produced  the  desired 
turbulent  boundary  layer,  and  the  flow  was  approximately  symmetric  about  the  centerline. 
The  collected  data  were  then  used  to  calculate  the  correct  placement  of  the  upper  tunnel 
wall  to  achieve  a  uniform  freestream  velocity  throughout  the  test  section. 

Apparatus 

The  Virginia  Tech  Boundary  Layer  Wind  Tunnel  was  utilized  for  the  tests  performed.  The 
tunnel  test  section  is  24  feet  in  length  with  a  width  of  3  feet  and  a  variable  height.  The  test 
section  is  divided  into  three  eight-foot  sections  numbered  from  the  contraction  as  sections 
one,  two,  and  three.  The  first  section  has  an  accelerating  flow  up  to  five  feet  from  the 
contraction  and  approximately  a  constant  inviscid  core  velocity  downstream.  The  second 
section  has  a  constant  height,  but  the  third  section  has  a  linearly  increasing  height. 

A  custom  made  static  pressure  probe  with  a  tip  diameter  of  0.035  inches  and  a  Prandtl 
style  Dwyer  Pitot  probe  N28F  with  a  diameter  of  0.125  inches  were  used  for  the  pressure 
measurements.  The  two  probes  were  connected  to  an  electronic  manometer 
manufactured  by  Setra.  The  electronic  manometer  has  a  maximum  pressure  rating  of  5 
inches  of  water.  It  is  calibrated  for  one  volt  per  inch  of  water.  Voltages  were  read  directly 
from  a  Hewlett  Packard  3478A  multimeter  with  a  displayed  voltage  resolution  of  five 
significant  digits.  The  static  probe  was  mounted  in  a  threaded  barrel  and  gear  mechanism 
that  was  vertically  traversed  by  a  stepper  motor  that  was  powered  by  a  Hewlett  Packard 
6200B  power  supply.  One  revolution  of  the  threaded  barrel  was  equivalent  to  a  0.025  inch 
vertical  move. 

Experiments 

The  first  set  of  measurements  that  were  taken  was  at  1 1  spanwise  stations  located  21 .3 
feet  from  the  contraction  of  the  tunnel.  The  tunnel  height  at  this  x  location  was  11.375 
inches.  A  static  probe  was  mounted  so  that  the  stagnation  ports  were  at  the  same  x 
location  as  the  static  port  of  the  pitot  probe  but  at  a  y  location  that  was  in  the  freestream. 
The  wind  tunnel  was  turned  on  and  set  to  a  velocity  of  approximately  32  m/s.  The 
freestream  pressure  was  obtained  from  the  static  probe.  Once  the  freestream  pressure 
was  acquired,  the  pitot  probe  was  connected  to  the  electronic  manometer.  The  pitot  probe 
was  then  positioned  at  the  lower  wall  of  the  tunnel  and  vertically  traversed.  Thirty-one 
points  were  logarithmically  spaced  throughout  the  boundary  layer,  and  the  dynamic 
pressure  was  recorded  at  each  point.  From  these  pressure  readings,  a  velocity  profile 
was  calculated.  A  velocity  profile  was  taken  first  at  the  centerline  and  then  at  ±6,  ±9,  ±12, 
and  ±15  inches  on  either  side  of  the  centerline.  The  ±15  inch  profiles  were  located  3 
inches  from  the  tunnel  sidewall. 
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A  single  velocity  profile  was  taken  at  the  centerline  in  the  second  section  of  the  wind  tunnel 
at  a  x  location  of  9.5  feet  from  the  contraction.  The  upper  wall  height  was  10.375  inches  in 
section  2.  Again,  dynamic  pressure  readings  were  taken  at  31  points  throughout  the 
boundary  layer. 

The  final  tests  were  pitot  measurements  performed  at  15  streamwise  stations  along  the 
tunnel  centerline.  The  stations  were  spaced  at  approximately  12  inches  beginning  from  the 
second  tunnel  section  to  the  test  section  exit.  The  freestream  pressure  was  recorded  and 
the  velocity  was  calcu  lated  at  each  station . 

Results 

From  the  first  set  of  experiments,  the  dynamic  pressure  was  recorded  at  31  points 
throughout  a  boundary  layer  for  1 1  spanwise  locations.  Utilizing  a  FORTRAN  program,  the 
pressure  readings  for  each  z  location  were  nondimensionalized  using  the  freestream 
pressure  reading,  and  from  these  values,  the  nondimensional  velocity  values  at  each  point 
were  calculated.  The  velocity  profile  at  the  centerline  was  compared  to  a  previous  hot  wire 
velocity  profile  performed  in  the  same  tunnel  at  approximately  the  same  x  location  and 
under  the  same  conditions.  The  graph  in  figure  1  shows  the  comparison  between  the  two 
experiments. 


Figure  1 :  Comparison  of  Hot  Wire  and  Pitot  Probe  Measurements 
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Figures  2  and  3  compare  the  velocity  profiles  for  the  positive  z  locations  and  the  negative 
z  locations  respectively.  The  data  indicates  the  flow  is  approximately  symmetric  about  the 
centerline. 
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Figure  2:  Boundary  Layer  Profiles  for  Positive  Z  Stations 
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Figure  3:  Boundary  Layer  Profiles  for  Negative  Z  Stations 


For  each  of  the  z  locations,  the  freestream  velocity  and  viscosity  were  calculated  using  a 
Turbo  BASIC  program.  The  displacement  thickness  (8*)  and  the  momentum  thickness  (0) 
were  calculated  for  the  lower  boundary  layer  using  the  following  formulas. 


8*  for  the  upper  boundary  layer  was  approximated  using  the  following  relationships. 


s=y\±=.99 

ue 

Eauation  3 

s;  A 

i=A 

Eauation  4 

<5„  =. 99^(1.01(5,  H^J} 

Eauation  5 

II 

Equation  6 

Where  the  subscript  /  denotes  the  lower  boundary  layer,  the  subscript  u  denotes  the  upper 
boundary  layer,  h  is  the  tunnel  height,  hf  is  the  height  of  the  freestream,  and  y  is  the 
vertical  coordinate  measured  zero  from  the  lower  wall. 

The  shape  factor  (hf)  and  Ree  for  the  lower  boundary  layer  were  calculated  for  each  profile 
using  the  following  equations. 


Equation  7 


Re*  = 


U6_ 

v 


Equation  8 


Table  1  compares  all  these  calculated  values  for  the  11  locations.  Figure  4  is  a  graph  of 
the  lower  boundary  layer  parameters. 


Table  1:  Boundary  Layer  Parameters 


z  (in) 

Ue  (m/s) 

5,  (in) 

8,*  (in) 

8u*  (in) 

Ree 

H 

-15.0 

30.86 

1 .663E-05 

7.425 

0.563 

0.0758 

0.4674 

3.837 

0.291 

23194.48 

1.2041 

-12.0 

30.82 

1.658E-05 

3.960 

0.582 

0.1468 

0.4513 

5.817 

0.854 

22431.24 

1.2885 

-9.0 

30.84 

1 .662E-05 

4.455 

0.569 

0.1278 

0.4397 

3.342 

0.427 

21828.05 

1 .2947 

-6.0 

30.82 

1 .677E-05 

5.940 

0.631 

0.1063 

0.4843 

3.837 

0.408 

23805.89 

1 .3033 

0.0 

30.61 

1.660E-05 

5.445 

0.622 

0.1143 

0.4745 

4.332 

0.495 

22149.86 

1.3113 

6.0 

30.78 

1.671  E-05 

4.950 

0.633 

0.1278 

0.4828 

5.322 

0.680 

22586.14 

1.3103 

9.0 

30.86 

1 .673 E-05 

4.950 

0.611 

0.1235 

0.403 

3.837 

0.474 

23088.80 

1.5171 

12.0 

30.79 

1 .672 E-05 

3.960 

0.527 

0.1331 

0.774 

18850.07 

1.1263 

15.0 

30.82 

1 .675 E-05 

8.415 

1.159 

0.392 
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Figure  4:  Lower  Boundary  Layer  Parameters 


The  boundary  layer  in  the  second  section  of  the  wind  tunnel  was  studied  in  the  second 
experiment.  Again,  the  freestream  velocity  and  the  viscosity  were  calculated  using  the 
Turbo  BASIC  program.  The  displacement  thickness,  momentum  thickness,  and  the  shape 
factor  were  calculated  for  this  profile  as  well.  From  the  momentum  thickness,  freestream 
velocity,  and  the  viscosity,  Ree  was  calculated  to  be  approximately  1 9,000.  Figure  5  is  a 
graph  of  the  data  taken  at  the  centerline  in  section  2  and  3  of  the  tunnel  as  well  as  the  hot 
wire  measurements  and  the  data  from  the  Danish-Netherlands  Windtunnel  (DNW)1.  The 
DNW  data  is  for  a  Reynolds  number  of  approximately  20,000. 
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Figure  5:  Comparison  of  mean  boundary  layer  profiles  for  various  X  locations 


The  friction  velocity  (UT)  was  calculated  using  the  following  formula. 


Equation  9 


Were  c,  is  the  coefficient  of  friction  calculated  by  the  Ludwieg-Tillmann  equation. 

c  -0  246x10-°W^/Y°268 
cf  U./40X11)  ^  /V  J  Equation  10 

Capital  Delta  (A)  was  calculated  using 


A  =  S*  y 

i/Cf 


Equation  11 


Table  2  lists  the  values  for  Uref,  cf,  A,  and  UT  for  the  11  spanwise  locations.  Table  3  lists  the  same 
values  for  the  second  section  centerline  location,  the  third  section  centerline  location,  the  hot  wire  data, 
and  the  DNW  data. 


Table  2:  Additional  Boundary  Layer  Parameters 


z  (in) 

Uref  (m/s) 

cf 

Ut  (m/s) 

A  (in) 

-15.0 

31.93 

0.002539 

1.100 

15.796 

-12.0 

31.97 

0.002245 

1.033 

17.355 

-9.0 

31.92 

0.00224 

1.032 

17.011 

-6.0 

31.90 

0.002159 

1.013 

19.209 

0.0 

31.86 

0.002143 

1.002 

19.010 

6.0 

31.87 

0.002137 

1.006 

19.355 

9.0 

31.87 

0.001623 

0.879 

21.462 

12.0 

31.88 

0.002872 

1.167 

13.908 

15.0 

31.91 

0.001984 

0.971 

36.801 

Table  3:  Comparison  of  Data  Sets 


Data  Set  (in) 

Uref  (m/s) 

cf 

Ux  (m/s) 

A  (in) 

Section  2 

31.96 

0.002415 

1.0857 

10.72 

Section  3 

31.86 

0.002143 

1.002 

19.01 

Hot  Wire 

32.00 

0.002217 

0.995 

18.09 

DNW 

20 

0.002333 

0.6819 

22.89 

The  third  set  of  experiments  gave  the  freestream  pressure  readings  for  15  locations  on  the 
centerline  of  the  tunnel.  From  each  value,  the  velocity  of  the  flow  was  calculated  at  that  x 
location.  Figure  6  shows  the  streamwise  velocity  variation  down  the  tunnel  centerline. 


Figure  6:Streamwise  Velocity  Variation 


This  experiment  was  used  to  check  the  uniformity  of  the  velocity  throughout  the  tunnel  test 
section.  It  was  desired  to  have  the  same  velocity  in  the  third  section  as  the  velocity  in  the 
second  section.  Using  the  results  from  the  first  two  experiments,  a  new  tunnel  height  for 
the  third  section  was  calculated  using  the  following  formula. 


h 


iy _ 


Equation  12 


Where  h,  is  the  initial  height,  l/,  is  the  initial  velocity,  Unew  is  the  desired  velocity,  and  8U* 
and  8|*  are  the  upper  and  lower  displacement  thicknesses  respectively.  The  new  tunnel 
height  was  calculated  to  be  approximately  0.26  inches  lower  then  the  original  height. 


Conclusions 


After  comparing  the  velocity  profiles  for  the  1 1  spanwise  locations  in  the  third  section  of 
the  wind  tunnel,  the  flow  appears  to  be  symmetric  about  the  centerline.  Also,  the 
roughness  plate  placed  in  the  first  section  of  the  tunnel  produces  the  high  Reynolds 
numbers  needed  at  the  end  of  the  tunnel.  Lowering  the  upper  wall  in  the  third  section  of 
the  tunnel  will  increase  the  velocity  to  the  desired  uniformed  value. 

Reference 

1  AGARD-AG-335.  ‘Turbulent  Boundary  Layers  in  Subsonic  and  Supersonic  Flow.”  Edited  by 
Professor  William  S.  Saric,  p.16. 
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Pressure  Distribution  on  Flat  Plate  near 
Appendage 

Abstract 


Measurements  were  taken  to  document  the  mean  surface  pressure  distribution  in  a  wing/body 
junction  flow  at  Ree  =  Ue6/v  (Reynolds  number  based  on  the  momentum  thickness,  6  )  of 

23,200.  These  measurements  yielded  the  time  mean  static  pressure  distribution  on  the  flat  test 
wall  and  the  surface  of  the  wing,  and  the  data  set  complemented  the  mean  velocity  and 
turbulence  measurements  acquired  by  the  Laser  Doppler  Velocimetiy  technique.  These 
measurements  were  also  used  to  check  the  flow  symmetry  around  the  wing. 


Description  of  Apparatus 


The  tests  were  performed  in  the  Virginia  Tech  Boundary  Layer  Wind  Tunnel.  The  tunnel  test 
section  is  about  29  feet  in  length  with  a  width  of  3  feet  and  a  variable  height.  The  test  section 
is  divided  into  two  eight-foot  sections  numbered  from  the  contraction  as  sections  one  and  two, 
and  a  section  numbered  three  which  is  about  13  feet  long.  The  first  section  has  an  accelerating 
flow  up  to  five  feet  from  the  contraction  and  approximately  a  constant  inviscid  core  velocity 
downstream.  The  second  section  has  a  constant  height,  but  the  third  section  has  a  varying 
height  that  has  been  suitably  adjusted  to  get  a  zero  pressure  gradient  2-D  boundary  layer.  The 
tests  were  performed  in  this  section. 

The  wing  model  used  in  the  present  measurements  is  a  NACA  0020  tail  with  a  3:2  elliptical 
nose.  The  model  has  the  dimensions  of  9  inches  in  height,  a  chord  length  of  12  inches,  and  a 
maximum  thickness  of  2.824  inches.  The  wing  model  was  aligned  to  the  flow  by  considering 
the  pressure  taps  at  the  maximum  thickness  point  ( t ),  on  either  side  of  the  wing  ( x/t  =4)770, 
y/t -1.7263,  z/t  ~0.5).  Equalizing  the  pressure  at  these  locations  was  sufficient  for  the 
model  to  be  aligned  with  the  flow  (details  are  in  the  results  section).  To  reduce  the  blockage 
effect  generated  by  the  presence  of  the  wing,  0.25  inch  thick  plates  upstream  of  the  wing,  x  < 
11.875  inches,  are  placed  on  the  side  walls.  The  36.25  inch  spanwise  distance  of  the  tunnel 
side  walls  at  the  wing  changes  to  35.8125  upstream  and  downstream  with  the  insert  plates. 


To  measure  the  wall  static  pressure,  a  plexiglass  plate  with  230  pressure  taps  around  the 
contour  of  the  wing  and  a  small  metal  plate  with  139  pressure  taps  at  the  nose  region  are 
placed  on  the  tunnel  floor.  A  drawing  of  the  pressure  plate  is  shown  in  figure  1. 


Fig.l  Pressure  plate 

The  wing  model  with  24  pressure  taps  on  its  surface  was  mounted  on  the  plexiglass  plate,  and 
a  Scanivalve  system,  CTLR2P/S2-S6  Scanivlave  Corp.,  was  connected  to  all  taps  on  the 
tunnel  floor  and  on  the  nose  of  the  model  using  Scanivalve  diaphragms  and  Tygon  tubing.  The 
diaphragms  were  connected  sequentially  to  the  system.  The  Scanivalve  system  has  48  to  1 
multiplexing,  and  the  sampling  frequency  is  1000/sec.  The  data  was  sampled  for  2  seconds 
with  a  settling  time  of  30  seconds.  A  Setra  239  pressure  transducer  with  a  calibrated  range  of 
0-15  inches  of  water  was  used  to  measure  the  reference  dynamic  pressure.  Another  Setra  239 
pressure  transducer  with  a  calibration  range  of  -2.5  to  2.5  inches  of  water  was  used  to  measure 
the  static  pressure  at  the  pressure  taps.  An  inclined  manometer  with  a  resolution  of  0.01  inches 
was  used  to  verify  the  calibration  of  the  two  pressure  transducers.  A  data  acquisition  board 
(DT2801)  in  conjunction  with  a  PC-AT286  was  used  to  acquire  the  data  from  the  pressure 
transducers. 


Procedure 


A  Prandtl  type  Pitot-static  probe  mounted  about  58  inches  from  the  end  of  the  upstream 
contraction  at  the  test  section  entrance  is  connected  to  the  reference  pressure  transducer.  The 
static  port  of  the  pitot  probe  is  also  connected  to  the  low  end  of  the  measurement  pressure 
transducer.  The  output  of  the  two  pressure  transducers  is  read  by  the  data  acquisition  board. 
The  pressure  readings  are  reduced  by  the  data  acquisition  computer,  and  the  pressure 
coefficients  are  calculated  using. 


P  —  P 

1  o  ref  rref 


The  subscript  ref  denotes  values  at  the  reference  station.  Pref  is  the  static  pressure,  Poref  is 
the  stagnation  pressure,  and  P  is  the  pressure  read  at  each  location. 

The  data  set  (shown  in  Table  1)  obtained  was  made  when  the  tunnel  was  run  at  a  32  m/s 
nominal  reference  velocity.  This  run  took  approximately  200  minutes,  and  this  included  the 
time  required  for  the  diaphragms  to  be  changed.  The  uncertainty  in  the  measurement  of  Cp  is 
over  the  domain  is  Cp  ±  0.018  (20  to  1  odds). 

Results 


Figure  2  shows  a  contour  plot  of  the  Cp  values.  Also  plotted  on  the  graph  is  the  data  taken  by 

Oilmen  and  Simpson  (1995)  for  the  same  wing/body  configuration  at  a  Re*  «  6000.  The 
contours  are  similar  to  each  other. 

Figure  3  shows  the  comparison  of  the  measured  values  with  the  potential  flow  values 
calculated  in  Oilmen  and  Simpson  (1995).  The  graphs  in  figures  4a  and  4b  are  plots  of  the 
data  taken  on  the  surface  of  the  wing  at  two  Y/T  locations.  The  data  collected  from  the 
pressure  taps  on  the  wing  and  from  the  pressure  plate  are  listed  in  Table  1. 


Fig  3.  Measured  and  calculated  Cp  contours.  (Contours  inside  wing  shape  are  meaningless 


but  show  contour  values) 


Oilmen  and  Simpson  (1995)  carried  out  LDV  measurements  at  eight  locations  in  the 
wing/body  junction  flow  at  Re^  ~  6000  (see  Fig.  5).  For  the  present  case  of  Ree  ~  23,200 
flow,  the  pressure  gradient  vectors  for  these  eight  locations  are  calculated  from  the  pressure 
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measurements.  These  values  are  in  Table  2.  The  values  indicate  that  the  at  the  same  stations, 
the  magnitudes  of  the  pressure  gradients  in  both  jc  and  z  are  lower  in  the  high  Reynolds 


number  case  as  compared  to  the  lower  Reynolds  number  case,  indicating  that  the  higher 
Reynolds  number  flow  subjects  these  stations  to  a  milder  pressure  gradient  as  compared  to  the 
lower  Reynolds  number.  As  mentioned  earlier,  the  wing  model  was  aligned  to  the  flow  by 
considering  the  pressure  taps  at  the  maximum  thickness  point  (t),  on  either  side  of  the  wing 
(x/t -0.770,  y/t  -1.7263,  z/t  ~0.5).  Equalizing  the  pressure  at  these  locations  was  sufficient 
for  the  model  to  be  aligned  with  the  flow.  This  is  seen  in  Fig.  6  where  the  Cp  values  are 
plotted  versus  z/t  values  at  either  sides  of  the  wing  at  constant  x/t  values. 


Fin  4h:  Cn  values  on  winn  at  Y/T=1  7?f> 


1.0 


.5  2.0  2.5 


-1 .0  -0.5  0.0  0.5 

x/t 

measurement  locations.  Full  arrows  wall-stress  direction.  Empty 
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Conclusions 


Pressure  measurements  were  mapped  in  a  wing/body  junction  at  a  Reynolds  number  of 
23,200.  The  measurements  also  yielded  the  pressure  gradients  at  the  eight  locations  where 
LDV  measurements  have  already  been  carried  out.  The  values  indicate  that  at  the  higher 
Reynolds  number  flow  these  stations  are  subjected  to  a  milder  pressure  gradient  as  compared 
to  the  lower  Reynolds  number  flow.  The  data  obtained  indicated  that  the  flow  is  symmetric 
around  the  wing,  and  this  is  shown  in  figure  6.  It  should  be  noted  that  taps  on  the  wing  at 
( x/ 1  =0.770,  y/t~ 1 .7263,  z/ 1  =0.5)  were  used  for  the  alignment  of  the  wing,  and  this  proved 
sufficient  for  the  flow  to  be  symmetric  elsewhere. 
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Uncertainty  Analysis 

The  uncertainty  in  the  measurement  of  Cp  is  calculated  by  the  method  outlined  by  Kline  and 
McClintock  (1953).  As  mentioned  earlier,  Cp  is  defined  as 

c  p~p"' 

The  measurement  of  Cp  with  the  two  transducers  (called  the  measurement  transducer  and  the 
reference  transducer)  were  calculated  using  the  reduced  expression  for  Cp  which  is 

£-  —  P  ~  Eisi  -  )mg as 

"  Po«f-P,.f  (A/>Lf 

It  may  be  noted  that  (AP)meas  is  measured  using  the  measurement  transducer  (range  of  -2.5  to 
2.5  inches  of  water)  and  (AP)re/  is  measured  using  the  reference  transducer  (range  of  0  to  15 
inches  of  water).  The  uncertainty  in  the  measurement  of  Cp  would  be  due  to  the  uncertainty 

in  the  measurement  of  these  two  quantities.  For  the  purpose  of  calculation  the  range  of 
pressure  sensed  is  taken  to  be  2.5  inches  of  water. 


Let  b  [(AP)meai  ]  be  the  uncertainty  in  (A P)meas  and  b  [(AP)re/  J  be  the  uncertainty  in  (A P)ref  , 
then  the  uncertainty  in  Cp ,  which  is  bCp,  is  given  by 


Y 

J 


Y 

J 


Uncertainties  in  the  individual  transducers 
(i)  Measurement  transducer  (-2.5  to  2.5  inches  of  water) 
Non-linearity :  ±  0.1  %  FS  (lull  scale) 

Zero  error  :  ±  0.4  %  FS 
FS  error  :  ±  0.4  %  FS 


(ii) 


Therefore,  <5 [(AP)m<ll J  is  equal  to  ±  Vo.0012  +  0.0042  +  0.0042  =±  0.0058  FS 
<5(AF)„ 


That  is, 


(AF), 


r meas  •  , 

lo  - 


±  0.0058 


Reference  transducer  (0  to  15  inches  of  water) 
Non-linearity :  ±  0. 1  %  FS 
Zero  error  :  ±  0.4  %  FS 
FS  error  :  ±  0.4  %  FS 


Therefore,  b  [(AP)meo?  ]  is  equal  to  ±  Vo.0012  +  0.0042  +  0.0042  «  +  0.0058  FS 
b  (A P). 

For  our  range  of  pressures  ,  —  J  works  out  to  be  ±  0.0058  x  15/5  ~  ±0.0174 


(A F), 


ref 


Hence,  the  net  uncertainty  in  C„  is  ±Vo.00582  +0.01742  ~±  0.018  which  is  ±  0.018 


Table  1.  DataSet 


Flow  Conditions 


Flow  temperature  (°C)  =  24.2 

density  (Kg/m3)  =1 .093452 

viscosity  (m2/s )  =  1.676351E-05 

Atmospheric  pressure  (Pa)  =  93300 

Velocity  of  undisturbed  free  stream  ( U  rej  in  m/s)  =31 .9035 1 

Reynolds  number  based  on  Uref  and  chord  of  wing  ( c )  =  580080.8 


Table  1  (Continued) :  Data  collected  from  the  wing 


ON  THE  WING 

X/T 

Y/T 

Z/T 

Cp 

0 

1.4607 

0 

1.007089 

0.28463 

1 .4607 

-0.3921 

-0.547441 

0.76991 

1.4607 

-0.5001 1 

-0.68712 

1.8357 

1.4607 

-0.4319 

-0.247373 

0.28526 

1.4607 

0.39244 

-0.465812 

0.77881 

1.4607 

0.50008 

-0.672856 

1.1327 

1.4607 

0.49014 

-0.429646 

1.4818 

1.4607 

0.4666 

-0.327168 

1.83 

1.4607 

0.43253 

-0.252858 

2.1769 

1.4607 

0.39014 

-0.193071 

2.5331 

1.4607 

0.33938 

-0.136104 

2.8771 

1.4607 

0.28434 

-0.084942 

0 

1.7263 

0 

1.007089 

0.28213 

1.7263 

-0.39078 

-0.542879 

0.76991 

1.7263 

-0.5001 1 

-0.696738 

1.1229 

1.7263 

-0.49063 

-0.42889 

1.4805 

1.7263 

-0.46671 

-0.32621 1 

1.8307 

1.7263 

-0.43245 

-0.249586 

2.1822 

1.7263 

-0.38944 

-0.190772 

2.5373 

1.7263 

-0.33874 

-0.122926 

2.889 

1 .7263 

-0.28234 

-0.086234 

0.27589 

1.7263 

0.38742 

-0.499354 

0.77102 

1.7263 

0.5001 1 

-0.694919 

1.8441 

1.7263 

0.43096 

-0.252828 

-0.31232 


-0.25 


-0.25  ] 


-0.25 


-0.25 


-0.25 


-0.25 


-0.25 


-0.25 


-0.25 


-0.18732 


-0.125 


-0.125 


-0.125 


-0.125 


-0.125 


-0.125 


-0.125 


-0.125 


-0.125 


-0.125 


-0.125 


-0.06267| 


0.125 


0.25 


0.375 


0.50035 


0.87571 


1.00041 


1.2504 


1.5004 


0 


0| 


0.062677 


0.125 


0.18732 


0.25 


0.31232 


0.375 


0.43732 


0.50035 


0.56303 


0.62571 


0| 


0.125 


0.18732 


0.25 


0.31232 


0.375 


0.43732 


0.50035 


0.56303| 


0.45136 


0.42881 


0.37924 


0.31585 


0.25518 


0.20877 


0.17175 


0.14368 


0.11407 


0.09317 


0.073581 


0.57041 


0.52595 


0.37816 


0.26008 


0.16081 


0.083339 


0.035148 


0.016238 


0.00173 


0.75035 


-0.33103 


-0.03001 
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Table  2.  Pressure  gradient  vectors  at  the  measurement  locations 


Station 

Ree  =23200 

Re5  =6000 

dCp 

w 

TC 

dCP 

TC 

*C, 

TC 

d^TJ 

TC 

i 

■EBH 

| 

2 

hemh 

■ 

3 

-0.2060 

0.2771 

| 

0.3312 

4 

0.1875 

ni'iiiiiB 

0.1454 

5 

-0.0587 

-0.4649 

-0.0649 

6 

-0.2643 

-0.3346 

-0.3609 

7 

-0.2544 

0.0122 

-0.4818 

8 

0.0395 

-0.1599 

0.0991 

-0.3161 

9 

0.0740 

-0.1094 

0.0926 

-0.1516 

